Realistic Exact Solution for the Exterior Field of a Rotating Neutron Star 
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A new six-parametric, axisymmetric and asymptotically flat exact solution of Einstein-Maxwell 
field equations having reflection symmetry is presented. It has arbitrary physical parameters of 
mass, angular momentum, mass-quadrupole moment, current octupole moment, electric charge 
and magnetic dipole, so it can represent the exterior field of a rotating, deformed, magnetized and 
charged object; some properties of the closed-form analytic solution such as its multipolar structure, 
electromagnetic fields and singularities are also presented. In the vacuum case, this analytic solution 
is matched to some numerical interior solutions representing neutron stars, calculated by Berti 
& Stergioulas P], imposing that the multipole moments be the same. As an independent test 
of accuracy of the solution to describe exterior fields of neutron stars, we present an extensive 
comparison of the radii of innermost stable circular orbits (ISCOs) obtained from Berti & Stergioulas 
numerical solutions, Kerr solution |3, Hartle & Thorne solution an analytic series expansion 
derived by Shibata & Sasaki 3 and, our exact solution. We found that radii of ISCOs from our 
solution fits better than others with realistic numerical interior solutions. 

PACS numbers: 04.20.Jb, 04.40.Nr, 95.30.Sf, 98.80.Jk, 02.30.Em 



I. INTRODUCTION 



Observed pulsars are assumed to be highly deformed ob- 
jects due to rotation. The spin frequencies of the 11 
nuclear-powered pulsars lie between 230 Hz (e.g. for 
the PSR 1845-19 pulsar)^ and 641-Hz (e.g. for the PSR 
B1937-^21 pulsar @, p), and within that range are 
marginally consistent with a uniform distribution (see 
0, 19 and references therein). The fast rotation is very 
important, not only in newly-born neutron stars which 
may undergo secular and dynamical instabilities (see Q 
and references therein) but also in strange stars. 

Studying Neutron Stars (NS) is interesting for several 
reasons. The quasi-stationary evolution of an isolated 
NS can be tracked considering equilibrium sequences in 
which the rest mass is constant while the angular mo- 
mentum varies [lol |. Such evolution may be driven by 
the adiabatic loss of energy and angular momentum via 
electromagnetic or gravitational radiation. On the other 
hand, the axisymmetric pulsations of rotating NS can be 
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excited in several scenarios, like core collapse, crust - and 
core - quakes or binary mergers, and could become de- 
tectable in either gravitational waves or high-energy ra- 
diation 111. Furthermore, the observational detection of 
such pulsations will yield valuable information about the 
equation of state of relativistic stars and therefore infor- 
mation about the properties and behavior of the matter 
at extreme conditions of high density {'-^ 10^^ gr/cm'^) 
and temperatures (around lO^K) jl2l |. In addition, new 
aspects of rotating NS have been revealed in about 1000 
pulsars. Some X-ray pulsars and some 7-ray pulsars have 
been detected in the past years. Among these new ob- 
jects, some exhibit quite different behaviour in their pulse 
periods fl^ . The measurement of the period and its time 
derivative yields evidence of ultra-magnetized stars, po- 
ssibly representing magnetars jl4j . 

In the last two decades, great advances have been made 
to understand the properties of astrophysical objects 
like White Dwarfs, Black Holes or NSs, in the frame 
of General Relativity. Particularly, from the numerical 
approach, several physical properties of NS such as en- 
ergy density, inertia moments, equatorial radii and oth- 
ers are known through the numerical solution of the full 
Einstein equations, assuming different equations of state 
for NS and applying algorithms, specially the Komatsu- 
Eriguchi-Hashisu self consistent field-method jlSi] and its 
variants (e.g. in 2005 by Ghezzi Gusakov, Yakovlev 
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& Gnedin [T^l and Berti, White, Maniopoulou & Bruni 
[|. In 2004 by Berti & Stergioulas Q and Yoshida & 
Eriguchi 13- In previous years by Stute & Camenzind 
in 2003 19\, Manko, Mielke & Sanabria-Gomez in 2000 
|20j. SibgatuUin & Sunyaev i21il in 1998, Bocquet Bonaz- 
zola, Gourgoulhon & Novak l22|. in 1995, and by Cook, 
Shapiro & Teukolsky in 1994ll3 and 1992 [H. and 
references therein). On the other hand, from the ana- 
lytic close approach, to construct axisymmetric station- 
ary exact solutions of the Einstein-Maxwell equations 
with physical parameters defined from the beginning is 
not a problem now, because of the development of the 
powerful integral method by SibgatuUin (see and 
|25|1 which allows to do this by considering the choosing 
of the Ernst's potentials on the symmetry axis. 

The aim of this paper is to present a new stationary ax- 
ially symmetric exact solution to the Einstein-Maxwell 
system as a model for the exterior field of a NS, not only 
for rapidly but also for slowly rotation rates. A consis- 
tent analytic closed representation of the exterior space- 
time around a rapidly rotating NS is desirable for many 
reasons, among them: (i) if we have an analytic closed 
form for the metric, the computation of the stationary 
properties of the space-time is less difficult (for example, 
geodesies in the exterior of NS could be studied analyti- 
cally [23 ; one could also find closed form expressions 
for the radii and frequencies of the innermost stable cir- 
cular orbits (ISCOs) ^ 0, etc.), (ii) it could be useful 
for studying the dynamical properties of the space-time, 
like gravitational wave emission and chaotic trajectories 
of particles around it. Furthermore, having an analytic 
solution could simplify the calculation of properties of 
accretion disks, epicyclic frequencies 21] and so on. 

Some analytic models for NS have been constructed and 
studied. Since the seminal works of Hartle [2^ and 
Hartle-Thorne '3], the analytic structure of the space- 
time outside a slowly rotating NS was generally associ- 
ated with the Kerr solution because of its simplicity and 
accuracy with numerical data obtained from numerical 
interior solutions. Nevertheless, Kerr solution does not 
present an accurate fit with the expected values of some 
of the physical properties of realistic stars in the regime of 
rapid rotation rates (like the radii of ISCOs, see |ll [l0l[l9| 
and the mass-quadrupole moment '2^). SibgatuUin & 
Sunyaev adjusted the free parameters of the exact 
solution by Manko, Martin, Ruiz, SibgatuUin & Zaripov 
in [s^l with the numerical data for the ISCO and the 
gravitational redshift of Cook, Shapiro & Teukolsky [Toj . 
They analyzed the effects of the mass-quadrupole mo- 
ment of a rapidly rotating NS on the energy release in 
the equatorial layer on the surface of the accreting star 
and in the accretion disk, using the normal sequences of 
Equations of State (EOS) A [sil, AU [p and EPS 



the charged, magnetized generalization with arbitrary 
quadripolar deformation of the Tomimatsu-Sato 6 = 2 
solution |34j written in rational functions, and claimed 
this solution as a possible model to describe the exte- 
rior field of a rotating NS. Then, Stute & Camenzind 
[1^ decided to study that solution to describe the grav- 
itational field of a rapidly rotating NS. They matched 
the solution by Manko et al. with the numerical inte- 
rior solutions from Cook et al ^3 by a similar procedure 
to the one used in concluding that the accuracy of 
that solution is high in the case of the redshift, but poor 
describing the radius of the marginally stable orbit (the 
matching was made using non invariants local properties 
of the solution, as it was pointed out in Qj). However, 
they considered that Manko et al. solution was a great 
improvement compared with Hartle & Thorne solution 
in the regime of fast rotation. 

Recently, the full Einstein equations were solved in a nu- 
merical approach by Berti & Stergioulas (B&S) 0, to 
determine the NS space-time along of sequences of con- 
stant rest mass for selected EOS (denoted as EOS A [sif . 
EOS AU [13, EOS EPS jH, EOS L jH and EOS APRb 
[sef). They matched the Manko et. al. solution "^^l to 
the numerical solutions imposing the condition that the 
mass-quadrupole moment of the numerical and analytic 
space-times be the same, concluding that the matching 
condition can be satisfied only for very rapidly rotat- 
ing stars. This affirmation was based on the fact that 
this exact solution does not reduce continuously to the 
Schwarzschild one when the rotation vanishes and, ac- 
cording to 1] , the mass-quadrupole moment in the non- 
rotating limit (angular moment equals to zero) is very 
large for this fimit, i.e., Q = — m(m^ + 6^)^/(4(m^ — 6^)) , 
where m is the mass and 6 is a parameter related with 
the arbitrary mass-quadrupole moment respectively. For 
that reason, it is not possible that the solution by Manko 
et al. can describe NSs in the regime of slow rotations. 
Besides, B&S suggested that for intermediate rotation 
rates could be used the exterior approximated analytic 
solution by Hartle & Thorne ^ , valid to second order in 
the rotation rate. 

Following their own suggestion, in 2005 Berti et al. 
compared the Hartle-Thorne Q slow-rotation approxi- 
mation keeping terms up to second order in the stellar 
angular velocity, and Manko et al. solution |20j | again 
with their numerical solutions to the full Einstein's equa- 
tions. Using the same matching procedure presented in 
0, they found that the Hartle-Thorne approximation 
offers very good predictions for the corotating i?+ and 
counterrotantig i?_ ISCO radii, with R± accurate better 
than 1% even for the fastest millisecond pulsars. At these 
rotational rates the accuracy on the mass-quadrupole 
moment is ~20%, and better for longer periods. 



In 2000, Manko, Mielke & Sanabria-Gomez ||2^ derived 
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served from pulsars motivates the study of the interior 
and exterior electromagnetic fields of rotating NS, we use 
the same procedure exposed in 0, i.e., the matching of 
the multipole moments in the vacuum case, because of 
the limited information available about electromagnetic 
properties in numerical interior solutions to Einstein- 
Maxwell system. However, considering that the study 
of electromagnetic fields in strongly curved space-times 
has been the subject of past and recent interest, we also 
present the closed form and some properties of the elec- 
tromagnetic field of our solution. 

This paper is presented in the following order. In sec- 
tion 2, we present the six-parametric solution including 
issues about its construction. In section 3, we summa- 
rize the properties of this solution through its multipolar 
structure, limiting cases, electromagnetic fields and sin- 
gularities. In section 4, we match our analytic solution 
with the interior numerical solutions by Berti & Ster- 
gioulas (hereafter B&S) and present comparisons of the 
radii of ISCOs of several numerical and exact solutions. 
Finally, concluding remarks about the paper are written. 



II. SIX-PARAMETER SOLUTION 



In order to construct a realistic exact solution to the 
Einstein-Maxwell equations, we must consider some 
physical assumptions. In Newtonian Theory of Gravi- 
tation all equilibrium states of isolated self-gravitating 
fluids (i.e. stellar models) must have reflection symme- 
try through a plane which is perpendicular to the rota- 
tion axis of the star [s^l; according to that idea, it has 
been conjectured that stationary general relativistic stel- 
lar models must have reflection symmetry as well |38) . 
For the stationary axisymmetric casg_a simple form of 
the metric was given by Papapetrou 39], that is 

ds^ = -f{dt ~ ujd(t>f + /"^ [e^^{dp^ + dz^) -f p^dct)"^] , 

a) 

where /, lo and 7 are functions of the quasi-cylindrical 
Weyl coordinates (p, z) . 

Using the above line element, the Einstein-Maxwell equa- 
tions can be reformulated, via Ernst's procedure 40], in 
terms of two complex potentials £ and $ as follow: 



{'R.eS + \^^)\/^£ = (V5 
(Re£:-Kl$l2)V2$ = {V£ 



2$*V$) • Vf, 
2$*V$)-V$. (2) 



In this paper, we solve the Ernst's equations (0) with 
the aid of the SibgatuUin's method il p], according 
to which the complex potentials £ and $ can be calcu- 
lated from specified axis data e(z) := £{z,p = 0) and 



<E>(z, p = 0), by the integrals 

1 e{Of^{<^)dcr 



£{z,p) 

$(z,p) = - / 
I" J-i 



TT 7-1 Vl - O-^ 

1 f' f{Ofi{a)da 



(3) 



The unknown function fJ.{cr) must satisfy the singular in- 
tegral equation 



' p{aM0+i{v) + 2f{0f{f^)]da 
-1 (cr — t)V1 — cr^ 

and the normalizing condition 
''^ fi{a)da 



= 



-1 



(4) 



(5) 



where ^ = z + ipa, rj — z + ipr, p and z being the Weyl- 
Papapetrou quasi-cylindrical coordinates, ct, r S [—1,1], 
e{rj) := e{ff), fijf) := f{fj) and the overbar stands for 
complex conjugation. 

With the purpose to construct reflectionally symmetric 
solutions, the Ernst's potentials on the symmetry axis 
e(z) and f{z) must satisfy the following conditions 

e{z)e{-z) = l and /(z) = -/(-z)e(z) , (6) 

when the electromagnetic potentials are even functions 
of z One choice of Ernst's potentials satisfying the 
above considerations is l431 



e(z) 



z^ — z^(m + ia) — kz + is 
z^ + z'^{m — ia) — kz + is ' 

v2 



qz- 



ipz 



z'^ + z'^im — ia) — kz + is 



(7) 



where m denotes the gravitational mass of the source, a 
its specific angular momentum, q its electric charge, and 
fc, s and p parameters related with the mass-quadrupole 
moment, the current octupole and the dipolar magnetic 

We have in- 



moment respectively (see Subsection IIII A|l 
eluded the parameter s following the suggestion given 
in 1!) regarding to the Manko et al. solution, whose 
accuracy is not good for the radius of the marginally- 
stable-orbit, which is more sensitive to higher multipole 
moments of the space-time. It seems to show that higher 
multipole moments ought to be included, especially for 
rapid rotation and stiffer equation of state, to improve 
the freedom to fit the exterior gravitational field of neu- 
tron stars with numerical interior solutions. Moreover, 
this parameter could aid us to understand the influence 
of the differential rotation in the dynamics properties of 
the surrounding space-time. 
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The potentials jT)) can be written in an alternative way, with 
we mean 



1=3 



i=3 



with 



(/3,-/3,)(/3,-/30 ' 



/. = (-1) 



Then, using and Q , following a similar procedure to 
the one used in ^^li we obtain the Ernst potentials and 
the metric functions in whole space-time: 



£ = 



A + B 
A-B' 



$ = 



C 



A-B' 



(8) 



^ij k ^ 
Cij = 

+ 



9^j 



Pij k 



\ ^) ■'■Hjk l\mn i 

( 1) -^l\mnp 1 

^ n\pl 

^3 ^12|/mn) i 

( 1 ) ^ij {^l ^m\np ~ '^m -^^7,1^/ 

'^p\lm ^ ^p -^/|mn) ; 
i-iyDi[f{ai)Hm\nps - /("m) -ffn|psi 

63 i 12\lmnp) 7 
63 ^' 12 limn) I 



and 



/ = 



- SB + CC 



=27 _ 



- BB + CC 



lui[{A + B)H -{A + B)G - CI] 
~ AA-BB + CC ' 

where 

A = ^ aijkrirjVk, B= ^ bijrirj, 

l<i<j<k<e l<i<j<6 

C ^ J2 ^^J'^' ' H = zA-{fii+l32+ P3)B 
i<i<j<& 

+ ^ kn rj rk+ ^ {ai + a/) bij n rj , 

l<i<j<k<e l<j<j<6 

G = -{pi+/32 + f33)A + zB+ 9ionrj 

l<i<j<6 

+ ^ {ai + aj+ak)aijkrirjrk, 

l<i<j<k<6 

I = (/i + /2 + f3)iA -B) + + /52 + /53 - z) C 

6 

+ X] + y^K^'i 

l<i<j<fe<6 i=l 
l<i<j<6 



Ay- 
Ay fe 



1^ 



(ai - Uj) Di Dj , 

(a,; - aj){ai - ak)iaj - at) Di Dj Dk 
1 

(a; - (3i){ai - P2){ai - /Ss) ' 

12 1 mn {am) Ai2|„Z 

-^3(0:™) Ai2|im , 

Hi{an) H,n{ap) - Hi{ap) i?„i(a„) , 

fl* fk 



3 

((3* - Pk){an - Pk) ' 



^ |mnp 



Zm|nps 
- /m|nprs 



Hs{ai) Si2\m7ip , 

^l\mnp ^~ ^m\npl ^~ ^7i\plm ~^ ^p\lmn i 
f{an) ^hn\ps + /(O^p) A/^|sn + /(«s) A/^|^p , 

rsn 



being a's the roots of the SibgatuUin's equation [2f 



eiz) + e{z) + 2f{z)f{z) = 0. 



(9) 
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III. PROPERTIES OF THE SOLUTION 



A. Multipolar structure 



This solution describes the exterior gravitational field of a 
reflectionally symmetric source. Their first four relativis- 
tic mass multipolar moments and first two electromag- 
netic multipolar moments are arbitrary. The multipo- 
lar moments are calculated using the Fodor-Hoenselaers- 
Perjes procedure (43i | following the parametrization given 
in 0, which has the advantange that the real parame- 
ters m corresponds to the total mass M, a to the total 
angular moment per unit mass (a — J/m, being J the to- 
tal angular moment) and q to the total electric charge Q; 
while fc, s and represent the mass-quadrupole moment 
M2, current octupole 5*3 and magnetic dipole B through 
the following formulas: 

A/2 — mk — ma^ , 
S3 — 2imak — iam'^ — ims, 
B = ifi + iaq. (10) 

The solution also has electric charge quadrupole moment 
Q2, which is not independent because it depends on the 
other multipoles, given by the formula 

Q2 = —afi — a'^q + kq , (11) 

and permits that the solution has electric field even when 
q vanish (induced electric field by rotation [isf'). 

For NSs the mass-quadrupole moment is, sur pris ingly, 
approximated by the simple quadratic relation |29l | 

M2^-ciM,E0S)^, (12) 

where the constant c — c{M, EOS) depends only on the 
mass M and the equation of state (EOS) for the interior 
of the NS. So, the formulas Q are successfuU to fit M2 
of the exact solucion with data of numerical solutions 
because of the form of M2 and the freedom to set the 
numerical value of the real parameter k. In addition, 
the absence of the gravitomagnetic parameter (monopole 
of angular momentum) leads to achieve the asymptotic 
flatness condition of the solution, necessary to describe a 
compact object. 



B. Limiting Cases 



A special feature of this solution is that in rigid rotation, 
i.e., with s = 0, it contains the following well-known 
limiting cases: 



1. The stationary vacuum case g = 0, /z = with non- 
vanishing arbitrary deformation parameter k = 
— l/4(m^ — a^) represents the Tomimatsu-Sato S — 
2 solution |34| with the mass- quadrupole moment 
M2 = -l/4(m3 - j2/m). 

2. The magnetostatic limit q = 0, a — represents the 
massive magnetic dipole solution of Bonnor . 

3. The Kerr-Newmann space-time is obtained with 
fi^k^O. 

4. The stationary vacuum case g = 0, = with non- 
arbitrary deformation, i.e., with k — represents 
the Kerr space-time. 

It is important to emphasize that the present six para- 
metric solution generalizes the Kerr and Kerr-Newmann 
solutions, which are the most studied stationary axisym- 
metric solutions. 



C. Electromagnetic Fields 



The strong magnetic field is one of the most important 
features of NSs and pulsars. The last ones have non- 
aligned dipole moment with the rotation axis causing loss 
of energy and angular moment, so the star radiates elec- 
tromagnetic and gravitational waves away and decrease 
its velocity of rotation, configuring a clear example of 
a nonstationary and nonaxisymmetric system. On the 
other hand, the stationary axisymmetric case is charac- 
terized by the fact that the magnetic dipole moment is 
aligned with the axis rotation of the star and the mag- 
netic field is poloidal (i.e. the magnetic field lying in the 
meridional planes p^lV Our exact solution is stationary 
and axisymmetric and, therefore, describes this kind of 
NSs. 

For the description of the electromagnetic properties of 
the solution, we must to calculate the electric and mag- 
netic fields produced in the surrounding space-time. The 
reason to introduce electric field in the solution, in spite 
of the fact of global neutrality of the astrophysical ob- 
jects, is that an electric field is induced by the rotation 
of the magnetized objects, as Eq. Hll() shows. To have a 
closed form the electric and magnetic fields is important, 
for example, in order to study the accretion of charged 
particles around of the NS, to study geodesic motion of 
test particles around a NS and so on. The electric and 
magnetic fields can be calculated using the expressions 

= F^pu^ , = -le^fj^'^F^su^ , (13) 

where Fap is the electromagnetic field tensor Fap = 
2 A[i3.a]i = (0, 0, A^, —At) is the electromagnetic four- 
potential, eap-yS is the totally antisymmetric tensor the 
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positive orientation with norm eap-^sf-"'^^^ = —24 
and Ua is a time-like vector (in the case of fluid, this 
vector is the four- velocity of the fluid). For a congruence 
of observers at rest in the frame of (Q), the four-velocity 
is defined by the time-like vector 



(1/77,0,0,0), 



(14) 



so the vectorial fields have components in the p and (p 
directions. The electric field is given by the following 
expressions: 



En 



and the magnetic field by 

f3/2 



B„ 



pe 



27 



p/2 

^—^{~u;Af,p + A^,p). 

pc 



(15) 



(16) 



(17) 



The At potential is the real part of the electromagnetic 
Ernst potential $ given by , and the potential A^ can 
be calculated as the real part of the Kinnersley potential 



IC = A, 



+ i 



iA'^ 48], which can be obtained using the 



SibgatuUin method and can be written as 



IC 



/(/1 + /2) 
A-B 



(18) 



Thus, the closed-form expressions for the electric and 
magnetic fields are 



En 



E. 
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-1 



\A-BW\A\^-\B\^ 

|2T-r6 



\c\' 



:Re< 



\A-BW\A\^-\B\^ + \C\^ 



C.. 



Cln(A-B), 



A-B 

C,-Cln(A-g),, 
A-B 



(19a) 
(19b) 



B„ 



Im(A - B)H +[A- B)G - Cl]^ , \K\^ ^\A\^ - \B\^ + \C\^ ULi \ (/i + h) [/. - IMA - B),,] 



p\A-B\ 



=- 



p\A-B[ 



lm{A-B)H+{A-B)G-CI]^ , \K\^ V\W^W+W ULi 



p\A-B[ 



-E^ 



p\A-B[ 



■ Im< 



A-B 



(20a) 



(/1+/2) [I,p-IHA-B)^p\ 



A-B 



(20b) 



r 



In order to show that our solution have the appropriate 
features of realistic NSs, we have to use typical values 
of NS for their physical parameters, like mass, rotation, 
deformation, charge and magnetic dipole moment. NS 
are compact objects that possess strong magnetic fields 
{B ~ lO^^j-)^ nevertheless, in our solution the parameter 
related with the magnetic field is the magnetic dipole mo- 
ment p. In pb] Manko et. al. showed how the magnetic 
dipole moment of their exterior solution depends on the 
values of the parameters of the interior solution, such as 
the matter density, the current function and the stellar 
radius. They obtained the above result by comparison 
between their magnetic potential A^ at infinite and the 
Ferraro's potential this became a test for the weak 



field because the Ferraro's solution is an approximate one 
that have into account deformation only at first order 
around the spherical symmetry According to this, 
we must to use typical values given for realistic interior 
solutions. Bocquet aZ. 22] calculated interior solutions 
for the Einstein-Maxwell system using the typical value 
of IQ^^Arn^ for the magnetic dipolar moment of the NS 
in S.I. units. In natural units {c = G = po = = 1), 
which have been used in our paper, the magnetic dipole 
moment has a order of 



l^natural ^ o f-^S.I. i 



(21] 



resulting a typical value for the p parameter of our solu- 
tion of order lO^Km^. 
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In the Fig. n]we have plotted the force hnes of the mag- 
netic field and the isopotential lines of the induced elec- 
tric field (At =constant) for one possible model of NS 
{Mq = lA7663Km). The increasing of the separation 
between consecutive force lines indicates, as is usual, that 
the electric field and the magnetic one decreases while 
the distance increases. As we expected, in absence of 
monopole electric charge distribution, the electric field 
is quadrupolar. In addition, the same figure shows that 
the axistationary properties of the magnetic field appear: 
aligned magnetic dipole moment and dipolar structure 
(poloidal field). 



D. Singularities 

The singularities are solutions of the equation 

A-B = 0, (22) 

i.e., denominator of the Ernst potentials given by 
equals to zero. In the Fig. |21we have plotted the typical 
shapes of stationary limit surfaces (/ = 0) for specific 
sets of realistic parameters given in Ref. Q. The lo- 
cation of ring singularity appears with dots. The ring 
singularity lie on the stationary limit surface, which is a 
characteristic of the stationary vacuum solutions. 

For all stellar models there are no ring singularities out- 
side of infinite red shift surfaces, therefore the solution 
can represent the exterior field of compact objects. In 
addition, we show in this section, with an approximate 
calculation via Komar masses, that the ring singularities 
of solution are massless. The geometrical structure of the 
exact solution implies that one could consider it as super- 
position of subextreme and hyperextreme objects. The 
Komar mass Ali of a subextreme part can be calculated 
using the Tomimatsu formula |^ 



(a) 



(23) 



where Vt is the imaginary part of the Ernst potential e, 
and is the constant value of the metric coefficient lo 
on the segment [a<, q;>] of the z-axis. 
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FIG. 1: (a) Magnetic field force lines and (b) isopotential lines 
of the induced electric field for m = 1.742M0, a = l.OQQKm, 
k = -O.SSeifm^ s = 0.212J^m^ q = OKm, b = lOKm'^ . 



In the case of hyperextreme part, the Komar mass can- 
not be evaluated using the formula if^ . and one more 
general integral expression should be used (sj 




[pO-T^f),p-(^^,z]p=Podz+ /[p(ln/) 



n,pl^,^dp- J [p(ln/)„~c^f],p],=,,dpl(^4) 



where the integration is carried out over the surface of a 
cylinder enclosing the cut which joins the points a; and 
af, Zu and z; denote locations on the symmetry axis of 
the centers of upper and lower bases, po is the radius of 
the bases. We will illustrate the procedure followed us- 
ing the set of parameters showed in the Fig. |3 First, 
we shall calculate the mass of the NS via Komar masses 
and then we will compare the obtained result with the 
mass monopole to, to deduce that the mass of the ring 
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0.5 1 1.5 2 2.5 



FIG. 2: Stationary limit surfaces for m = 2.388M0 = 
3.526A"m, a = 2.264Km, k = -l.SOSKm^, s = 2.498Km^ . 



singularity is equal to zero. For the set of parameters 
showed in Fig. [3 the roots of the SibgatuUin equation 
© are m = 0.792, a2 = 2.245,^3 = -0.792, a4 = 
— 2.245, = — 1.406i,a6 = as, producing four subex- 
treme parts and one hyperextreme part. Using the for- 
mula H23|l . the mass of the subextreme contribution is 
nis = 1.523Km and using the formula H24|) . the mass of 
the hyperextreme contribution is to^ = 2.003/^771, then 
nis + rrih — S.526Km — m. We made the above calcula- 
tions for all data for all EOS used from Ref. 1] and the 
results were allways the same. 

As can be seen, the calculations were made in a numer- 
ical way for the different stellar modells; it brings us to 
give only approximate values for the mass of ring singu- 
larity. For that reason, the ring singularity arise with a 
mass equals to zero with precision up to 10~*, with no 
dependence of the equation of state used. 



IV. TEST OF ACCURACY OF THE 
SIX-PARAMETER ANALYTIC MODEL 



In order to demonstrate that our solution can represent 
the exterior gravitational field of a realistic NS, we have 
to do comparisons with results obtained from numeri- 
cal interior solutions and with others possible analytic 
models. We have chosen the Kerr solution, the expres- 



sions derived by Shibata & Sasaki (S&S) in Ref.*^ and 
Hartle- Thome's solution (H&T) truncated to second or- 
der of the expansion parameter e — Q/fl* (see Ref.p^ 
and 0). For the comparison with numerical models, we 
have chosen the solutions by B&S Rcf. IJ; the comparison 
procedure was taken from there, that is, we compared an 
invariant quantity of the space-time, the multipolar mo- 
ments, and other quantity which can be depending on the 
coordinates. The depending coordinates quantity could 
be the gravitational red shift or the radii of the Inner- 
most Stable Circular Orbits (ISCOs). For this work we 
have chosen to compare the radii of the ISCO because 
it is related to several astrophysical properties of rapidly 
rotating NSs, e.g. the accretion disk cannot be longer 
that radii of the ISCO, and this sets an upper limit to 
the Keplerian frequency of particles orbiting a star. 

We matched the exact solution with the numerical inte- 
rior solutions imposing the condition that the first four 
mass and current multipole moments, i.e. mass, angu- 
lar moment, mass-quadrupole and current octupole, of 
the exact solution and the numerical ones have the same 
value. In B&S, sources possessing electric charge or mag- 
netic dipolar moment were not considered, so in this work 
the comparison of the solutions was developed in the vac- 
uum case. 

The numerical solutions were calculated for five interior 
numerical solutions for the gravitational field of a NS, 
each numerical solution is described by a different Equa- 
tion Of State (EOS). The EOS used were denoted as A 
01, AU 0, EPS 33], L 35], APRb ^ (see B&S [l| 
and Cook et al. [ifl,] for details). These' numerical so- 
lutions were calculated for three different sequences of 
constant rest mass: 



• First sequence. The sequence corresponding to a 
canonical neutron star having gravitational mass 
M = 1.4M0 in the non rotating limit. 

• Second sequence. The sequence terminating at the 
maximum-mass model in the non rotating limit 
(maximum-mass normal sequence). 

• Third sequence. A supramassive sequence, i.e., 
a sequence which does not terminate at a non- 
rotating model. 



In order to distinguish when a neutron star is rapidly ro- 
tating or not, we shall use the dimensionless parameter 
j, defined b y ?' = J/M^, as it was used by SibgatuUin 
& Sunyaev Stute & Camenzind [13 and B&S Q. 
In spite of j is not a direct quantitative measure of the 
magnitude of the rotation rate of the stars, it gives us a 
measure of the angular momentum which relates propor- 
tionally to the angular velocity. 
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To clarify the meaning of the value of the j parameter lets 
consider an example. For the recently discovered pulsar 
which frequency is 716-Hz, the value of j depends on 
the EOS and the sequence of mass considered. Using the 
values reported in B&S for stars with similar frequency, 
we find that j value must be around 0.42 < j < 0.68. For 
example, for the EOS EPS in the first sequence of mass 
and assuming m ~ 1.430A/q and J ^ 2.961 Km^, the 
value of j must be very close to 0.67, and for the second 
sequence assuming m ~ 1. 831 Mq and J ~ 3.185 if m^, 
the value of j must be very close to 0.43. Although, it 
is clear that we can not present a reliable classification 
for rapidly rotating neutron stars based on j, we shall 
assume that a NS is considered rapidly rotating when 
j > 0.35. 



Innermost stable circular orbits 



It is well known that not all orbits around relativistic 
stars are stable. Eor non-rotating stars, the ISCO is lo- 
cated at a circumferential radius of Risco — 6A/ (see 
Ref. 1^ for a complete treatment of circular orbits in 
stationary axisymmetric spacetimes). Depending on the 
EOS and the mass of the star, the ISCO can be located 
outside the stellar surface. The rotation introduces a pre- 
ferred direction in the coordinate, so ISCOs around a 
rotating star belong to two families: the co-rotating and 
the counter-rotating one. 

A circular orbit in the equatorial plane is one for which 
p = const. The equation for geodesic motion along the 
radial coordinate p reads 



dp 



9t4, - 9tt 94,4, 



for co-rotating orbits (-R-i-), where ' indicates a partial 
derivative with respect to p. For the case of Hartle- 
Thorne's solution, the radii of the ISCO for co-rotating 
orbits (i?+) was derived by Abramowicz et al. in p^: 



R 



Hk.T 
ISCO 



— 6m 




/251647 , 3\ 2 

27^ ' 2401^^2,)^ 



/9325 3 , 



(27) 



where g is a parameter related with the cuadripolar de- 
formation and defined as q = Q/nv", where Q is the cu- 
dripolar deformation multipole. 

S&S derived a general representation of axisymmetric 
vacuum solutions (in the form of a series expansion of the 
physical multipole moments of the space-time) and found 
some approximate analytic formulas for the location of 
the inner stable circular orbit, the angular moment, and 
energy of a particle around a relativistic source In 
general, their formulas depend on the mass, angular mo- 
mentum, mass- quadrupole, current octupole and mass 
2^-pole moments and on upper order multipoles in the 
rotation parameter. Including all terms up to order 0(4) 
in the rotation parameter, they found the following equa- 
tion for the circumferential radius of the corotating ISCO 
(i?+): 



pS&S 
^ISCO 



(25) 



6to(1 - 0.54433j - 0.22619j^ + 0.17989g2 

- 0.23002j3 -f 0.26296jg2 - 0.05317^3 

- 0.29693/ + 0.44546j2g2 - 0.06249g^ 

-I- 0.01544^4 -0.11310j(73). (28) 



where E and L are the conserved energy and angular 
momentum per unit mass, determined by the conditions 
V = dV/dp = 0. Geodesies become marginally stable 
when d'^V/dp^ = 0, so it is 

c.' u."fp {2f - f'p) + J^n2f + (-r + f"f)p'] 
+ Lo'f^u'^f^ + f'p{2f - !'p)[2f + 2/'2p2 



- ,/p(4/' + f'p) + p{2f - f'p){3f'f 

~ ^f"f p + f"p' + f[f"p 

- c."/v/^'2/^ + /'p(2/-/'p)]} =0, 



In the previous expression, dimensionless parameters 
q2 — —M2/m^, (73 = —S^/ra^ and — Mi/m^ have 
been introduced. S&S adopted the approximation, q^ — 
0452 ^ because the numerical calculation of the 2^-pole 
term is very hard and difficult. The value of a4 lies be- 
tween and 2, we adopted a4 = 1 as it is generally used 
0. In the case of the Kerr metric, the approximate ex- 
pression for the location of the co-rotating ISCO up to 
order 0{i'^) is 



(26) 



R 



Kerr 
ISCO 



= 6M{1 - 0.54433j - 0.04630j^ 
- 0.02016^-0.01110^*). 



(29) 



In all cases, the location of the counter-rotating ISCO 
{R-) is obtained simply using the above equation and 
changing the sign of all the star's rotation multipoles. 
The obtained results from the direct comparison between 
the ISCOs are presented in the Fig. |3| 
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FIG. 3: ISCO for the EOS A sequence with constant rest mass corresponding to a maximum-mass model in the non rotating 
limit of 1.948Af0 for (a.l) the counter-rotating orbits and (a. 2) the co-rotating orbits. ISCO for the EOS AU sequence with 
constant rest mass corresponding to a nonrotating model of 1.578M0 for (b.l) the counter rotating orbits and (b.2) the co- 
rotating orbits. ISCO for EOS EPS sequence with constant rest mass corresponding to a ma^ximum-mass model in the non 
rotating limit of 2.226Mq for (c.l) the counter-rotating orbits and (c.2) the co-rotating orbits. ISCO for EOS L sequence 
with constant rest mass corresponding to the maximum-mass model in the non rotating limit 3.232Mq for (d.l) the counter- 
rotating orbits and (d.2) the co-rotating orbits. ISCO for EOS APRb sequence with constant rest mass corresponding to the 
maximum-mass model in the non rotating limit 2.672Mq for (e.l) the counter-rotating o rbits and (e.2) the co-rotating orbits. 



We can see that our solution presents a better adjust 
for the radii of the ISCO in comparison with the stan- 
dard analytic models used in the literature to describe 
the exterior gravitational field of a rotating neutron star, 
e.g. Kerr solution, Hartle & Thorne' solution (H & T' 



solution) and the Shibata & Sasaki' expressions (S & S' 
expression). In average, for all cases the accuracy of R+ 
is better than 2%. and for i?_ better than 1%. The fact 
that the six parametric solution fits quite well from the 
non-rotating limit to the rapidly rotations rates, in all 
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sequences of mass in the five studied equation of state, 
implies that it could be an appropriated analytic closed 
form model for the exterior gravitational field of a neu- 
tron star. 



V. CONCLUDING REMARKS 



We present a new stationary axisymmetric six-parameter 
closed-form analytic solution generalizing to Kerr- 
Newmann solution with arbitrary mass-quadrupole mo- 
ment, octupole current moment and electric charge. The 
analytic form of its multipolar structure and their electric 
and magnetic fields, are also presented. 

We calculated some properties of the solution, using re- 
alistic data of interior numerical solutions, in order to 
demonstrate that the solution could model the exterior 
field of a realistic Neutron Star. In the clcctrovacuum 
case, we used the magnetic dipolar moment given by Boc- 
quet et al. in ^22]; i^i the vacuum case, we matched it 
to highly-accurate numerical solutions by Berti & Ster- 
gioulas 1], imposing that the mass, angular momentum, 
mass-quadrupole and current octupole moment of the 
numerical and analytic space-times have the same value. 
We also showed that this six-parameter analytic closed- 
form model is an appropriated model for the exterior field 



of a slowly or rapidly Neutron Star based on comparisons 
of the radn of ISCO obtained with a) Berti & Stergioulas 
numerical solutions (2004) 0, b) Kerr solution (1963) 
, c) Hartle & Thorne solution "s"] , d) an analytic series 
expansion derived by Shibata & Sasaki (1998) \^ with 
our analytic model for all Equations of State given in 0| . 
However, the solution's accuracy should be tested also 
outside of the equatorial plane, and through the calcula- 
tion of other physical observables, which is a task for the 
future. 

The exact solution could be used in future studies of as- 
trophysical plasma, dynamic of geodesies and accretion 
of particles in the surrounding space-time of the Neutron 
Star. 
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